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TECHNICAL PAPER 


A NEW APPROACH TO STATE ESTIMATION IN DETERMINISTIC 
DIGITAL CONTROL SYSTEMS 

I. INTRODUCTION 


Books on modem digital control systems usually address the problem of controlling a continuous- 
time plant driven by a zero-order-hold with a sampled output, as shown in Figure 1. For example, see 
Reference 1 , page 126. A common solution to this problem is to estimate the state of the system at the 
sampling instant using a state observer and then feedback the estimated state [ 1 , p. 1 95] . However, the state 
observer has two undesirable characteristics. First of all, it is a dynamical system in itself and, hence, adds 
additional states, eigenvalues, and dynamics to the system. This can affect system stability. Secondly, as a 
consequence, the estimated state is normally an approximation to the true state and is usually not a good one 
early in the estimation process. This paper presents a new approach to state estimation which has neither of 
these problems. The new estimator adds no new states, eigenvalues, or dynamics to the system and if the 
plant parameters are known exactly, the estimated state is actually equal to the true state. Useful in the 
development of the new state estimator are some results to date for continuous-time plants driven by a 
zero-order-hold with sampled outputs. These are reviewed in Section II, prior to the development of the 
new state estimator presented in Section III. Section IV presents an example which illustrates the procedure 
to follow to completely design the new estimator. Section V contains the conclusions and final comments. 
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Figure 1 . Continuous-time plant driven by a zero-order-hold 
with instantaneous measurements. 


II. PRELIMINARY 


For the plant in Figure 1 , x.(t) is an nxl state vector, u(k) is an rxl control input vector, yi(k) is an mxl 
output or measurement vector, F is an nxn system matrix, G is an nxr control matrix, and Q is an mxn 
output matrix. Since y ( (k) = Q x(k), where k is the usual shorthand notation for time kT, yj(k) represents 
an instantaneous measure of the system at the sampling instant kT. Hence, the plant in Figure 1 will be 
regarded as having instantaneous measurements for outputs. It is well known that this system can be 
modeled at the sampling instants by the discrete state equations [1, p. 126] 




x(k+ 1) = A x(k) + B u(k) 


( 1 ) 


y T (k) = Ci x(k) 


(2) 


where 


<f>(t) =^-'[(sI-F)-‘] , (3) 

A = 4>(T) , (4) 


and 


B = 


T 

[J <f>(X) dx] G 


(5) 


cj>(t) is the nxn state transition matrix. A and B are the nxn system matrix and the nxr control matrix, 
respectively, for the discrete state equations (1) and (2). 

A and B can be determined analytically using equations (3) to (5). An alternative approach, which is 
also quite suitable for numerical computation, is as follows [2] . cj>(t) and f ' <{>(\)d\ can be expressed in the 
form of matrix exponential series as 


<f>(t) 


2 ^ 

i =0 i! 


(6) 


and 


/ 4>(\)d\ = 2 ? 


■i fi + 1 


F t' 


i=o (i+1)! 


respectively. From equations (6) and (7), 


<K0 = I + F [/ 4>(X)d\] 


(7) 


( 8 ) 
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T 

where I is an nxn identity matrix. Hence, f 0 <}>(» d\ can be determined using equation (7) with t = T and 
this result substituted into equation (8) to get <f>(T). With these results, A and B can be found using 
equations (4) and (5). 

Now consider the plant in Figure 2, which is a generalization of the one in Figure 1 . In addition to 
the instantaneous measurement vector yj(kT), the plant in Figure 2 has the measurement vector y F (kT) 
generated as follows. First, the continuous-time output _z(t) is sampled every T/N seconds. Every N 
samples are multiplied by the weighting matrices H j; j = 0,1,...,N-1, and then summed to generate the 
output y F (kT), every T seconds. Functionally, this is equivalent to passing the discrete measurements 
generated every T/N seconds through a multi-input/multi-output moving average (MA) process with 
coefficient matrices Hj,j = 0,1,..., N-l [3], The output of the MA prefilter is sampled every T seconds to 
generate y F (kT). Then, y F (kT) has subtracted from it E_ u[(k-l)T], where E_ is a constant matrix, to pro- 
duce the modified MA-prefiltered measurement vector y F (kT). In Figure 2, C F is a pxn output matrix and 
z(t) is a pxl vector. The weighting matrices H j; j = 0,1,..., N-l are each qxp. Hence, y F (kT) is a qxl in- 
termediate output vector and y F (kT) a qxl output vector. Since u[(k- 1 )T] is an rxl delayed input vector, E_ is 
a qxr matrix. 

Previously, Polites [4] showed that when 


N-i j(T/N) 

E_ = L Hj C F J 4>(\)d\jG , 

j = 0 o 


the discrete state equations for the plant in Figure 2 become 


(9) 


x(k+ 1) = A x(k) + B u(k) 


( 10 ) 


y-r(k) 


'yi(k)" 


'c," 

_yp (k) _ 


D_ 


x(k) 


where D_ is a qxn matrix given by 


N-l 

D_ = S Hj C F <|>(-j -L) . 

j=0 N 

An alternative expression for equation (9) is 


( 11 ) 


( 12 ) 


E_ = H (3 


(13) 
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Figure 2. Continuous-time plant driven by a zero-order-hold with instantaneous and 
modified MA-prefiltered measurements. 


where 


H = [H 0 1 H, | 


! h n .,] 


(14) 


and 


3 = 


Cf 

C f 



4»(A) d\] G 


irwi 


G 



r -(N-l)(T/N) 

[J 4>(\) d\J G 


(15) 


In equation (14), H is a qx(Np) matrix and in equation (15), p is an (Np)xr matrix. An alternative expres- 
sion for equation (12) is 


D_ = Ha 


(16) 


4 





where 


a 


c F <M0) 

C F <J>(- X) 
N 


(17) 


C F <j>[-(N-l)-^] 


In equation (17), a is an (Np)xn matrix. 

E_ and D_ can be evaluated analytically using equations (3) and (13) to (17). An alternative 
approach, which can be either analytical or numerical, is as follows. Let t = -j(T/N) where j = 0,1,..., N- 
1 and use equation (7) to determine f Q jrT/N) 4>(\)d\, j = 0, 1 , . . . ,N- 1 . Substitute these results into equation 
(8) to get 4>[-j(T/N)] , j = 0,1,... ,N- 1 . At this point, E_ and D_ can be found using equations ( 1 3) to ( 1 7) . 


III. THE NEW STATE ESTIMATOR 


A general block diagram of the plant and the new state estimator is shown in Figure 3. Observe that 
this is Figure 2 with 


C, = 0 


(18) 


Thus, if E_ is given by equation (13), the discrete state equations for the complete system in Figure 3 are 


x(k+l) = A x(k) + B u (k) 

09) 

y F (k) = [Ha] x(k) . 

(20) 


This follows from equations (10), (11), (16), and (18). Recall that a is an (Np)xn matrix. If Nps*n, or 
equivalently N^(n/p), and a has maximal rank (i.e., rank n), then (a T a) is positive definite and 
nonsingular. In this case, H can be given by the pseudo inverse of a or [5] 


H = (a T a)-' a T 


( 21 ) 
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In general, H is aqx(Np) matrix. However, when it is given by equation (21), it is an nx(Np) matrix. This 
implies that q = n in this case. Consequently, yp(kT), which is a qxl vector in general, is an nxl vector in this 
case. From equation (19) to (21), the discrete state equations for the system in Figure 3 become 


_x(k + 1 ) = A 2 L(k) + B u(k) 

(22) 

>p(k) = x(k) . 

(23) 


Hence, the measurement vector y F (k) equals the true state vector _x(k). For H to be given by equation (21) 
and E_ by equation (13), the matrices F, G, and C F must be known. Given H, the coefficient matrices Hj, 
j = 0, 1 , . . . ,N- 1 in the MA prefilter are determined by partitioning H into N nxp submatrices. This follows 
from equation (14). 

At last, the new state estimation scheme, and the procedure for selecting the parameters in it, can be 
summarized as follows. A general block diagram of the new state estimator and the plant whose states are to 
be estimated is shown in Figure 3 . The MA-prefilter coefficient matrices Hj, j = 0, 1 ,. . . . ,N- 1 and the matrix 
E_ are to be calculated so the discrete state equations for the entire system in Figure 3 are given by equations 
(22) and (23). Consequently, the output vector y F (k)does more than estimate the state vector xfk); it equals 
it! Realize that forHj, j = 0,l ,...,N-1 and E_to be calculated so this is true requires that F, G, andC F in the 
plant be known precisely. Recall thatC F isapxn matrix. Choose the length of the MA pre filter, N, so thatN 
3= (n/p). Knowing N and F, J^ (T/N> c|)(\)dX. and <}>[-j(T/N)] , j = 0, 1 , . . . ,N- 1 , can be calculated using either 
of the methods described in Section II. Having these and knowing G and C F , a and (3 can be determined 
using equations (17) and (15), respectively. Test a to be sure that it has maximal rank (i.e., rank n). If so, 
let H be given by equation (21). The matrix E_can now be determined using equation (13). Partitioning H 
into N nxp submatrices as in equation ( 1 4) reveals the M A prefilter coefficient matrices Hj, j = 0, 1 , . . . ,N- 1 . 
The estimator is now completely defined. 


PLANT ESTIMATOR 



Figure 3 . General block diagram of the plant and the new state estimator. 
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IV. AN EXAMPLE 


Consider the double integrator plant driven by a zero-order-hold as shown in Figure 4. The con- 
tinuous-time output is sampled every T/N seconds and input into the new state estimator along with the 
control input u(kT). The parameters in the estimator are to be chosen so the discrete state equations for the 
complete system are given by equations (22) and (23). 

Manipulating the plant in Figure 4 into the format of Figure 3 yields 


[F] 


0 1 
0 0 


(24) 


[G] 



1 


(25) 


and 


[C F ] = [1 01 


(26) 


Using equations (24) and (25) and the formulas presented in Section II, 


<Mt) 


1 t 

0 1 


(27) 


PLANT ESTIMATOR 



Figure 4. Plant and new state estimator for the example. 
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(28) 


/ l 4>(X) 
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0 


t- 
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A 


1 

0 


T 

1 


(29) 


and 


B 


T" ‘ 

y 

T . 


(30) 


Since C F is a pxn matrix, it follows from equation (26) that p = 1 and n = 2. Hence, the requirement to 
select N so that N 5= (n/p) can be satisfied by letting 


N = 4 


(31) 


Using equations (15), (17), (25) to (28), and (31), a and (3 are found to be 


a 


1 

1 

1 

1 


0 

T 

4 

21 

4 

3T 

4 


and 


P = 


0 

j2 

32 

4T 2 

32 

9T 2 

32 


(32) 


(33) 
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respectively. In equation (32), it is apparent that the rank of a is maximal (i.e., rank = 2). Consequently, H 
can be given by equation (21). Substituting equation (32) into (21) yields 


H = [H 0 | H, j H 2 j H 3 ] = 


(■j i (“) I OiO 


(34) 


where the columns of H are the coefficient matrices for the MA prefilter by virtue of equation (14). Next, E. 
can be calculated using equations (13), (33), and (34) and is found to be 



The estimator is now completely defined. Using equations (32) and (34), one can verify that 


(35) 


Ha 


1 

0 


0 

1 


Hence, the discrete state equations for the entire system in Figure 4, with the estimator parameters defined 
by equations (31), (34), and (35), are given by equations (22) and (23) where A and B are defined by 
equations (29) and (30), respectively. 


V. CONCLUSIONS 


A new state estimator for deterministic digital control systems has been developed which offers two 
distinct advantages over the widely-used state observer. First of all, it contributes no additional states, 
eigenvalues, or dynamics to the system and, consequently, will not alter the stability of the system. 
Secondly, if the plant parameters are known exactly, the estimated state equals the true state and is not just 
an approximation to it. A disadvantage of this scheme is that the plant output must be sampled at a rate 
faster than the state observer’s. Fortunately, this is not the problem it would have been 20 years ago, con- 
sidering the speed of today’s digital computers. Still, if the sample rate on the output is so fast the computer 
solving the feedback control algorithms has a problem solving the estimation algorithms also, the following 
implementation can be used. A microprocessor can be dedicated to sampling the plant output and solving 
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the estimation algorithms to the point where the intermediate output ^F(kT) is generated. This can be sent, 
every T seconds, to the control computer solving the feedback control algorithms. There, the state estimate 
yp(kT) can be calculated by subtracting E_ u[(k-l)T] from y F (kT). This result can then be input into the 
feedback control algorithms to generate the control input u(kT). This process is repeated every T seconds. 

If the work in this paper is extended, the recommendation is to investigate the robustness of the new 
state estimator and see how it compares with the state observer’s. Specifically, the following questions 
should be addressed. What effect do modeling errors in the plant have on the new estimator and how does 
this compare with the state observer? What effect do plant process and measurement noise have on the new 
estimator and how does this compare with the state observer, or even the Kalman filter? How can the 
robustness of the new estimator be improved? Increasing the length of the MA prefilter, N, may be one 
possibility. Catenating the new estimator with a state observer or a Kalman filter could be another. This 
could produce a composite estimator which is better than either individual estimator by itself. 
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